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Abstract. Using simulations of isotropically forced helical turbulence the contributions to kinetic and magnetic alpha effects 
are computed. It is shown that for the parameter regimes considered in an earlier publication (Brandenburg & Subramanian 
2005, A&A 439, 835), the expressions for isotropic and anisotropic alpha effects give quantitatively similar results. Both 
kinetic and magnetic alpha effects are proportional to a relaxation time whose value, in units of the turnover time, is shown 
to be approximately unity and independent of the magnetic Reynolds number. 
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1. Introduction 

In a seminal paper, Pouquet al al. (1976) showed that in the 
nonlinear regime the alpha effect in mean field magnetohy- 
drodynamics is no longer governed by the kinetic helicity 
(Steenbeck et al. 1966), but there is an additional contribu- 
tion from the current helicity, so 
„, _ l, 



tT —U) ■ U 



Pa 1 3 



b) 
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where r is a correlation or relaxation time, u is the small scale 
velocity, u> = V x u is the small scale vorticity, b is the small 
magnetic field, and j = V x b//io is the small scale current 
density. Overbars denote some suitable form of averaging. 
Equation (T} has been used to explain catastrophic (magnetic 
Reynolds number dependent) quenching of the alpha effect 
in the nonlinear regime (Gruzinov & Diamond 1994; Bhat- 
tacharjee & Yuan 1995; Field & Blackman 2002). Techni- 
cally, the j ■ b term arises naturally when the t approximation 
is used (Kleeorin & Rogachevskii 1999; Radler et al. 2003; 
Blackman & Field 2002, 2003; see review by Brandenburg & 
Subramanian 2005a). 

In a recent paper, Brandenburg & Subramanian (2005b, 
hereafter BS05) presented results of numerical simulations 
that demonstrate the rise of the j ■ b term with magnetic 
Reynolds number in the presence of a finite imposed mag- 
netic field, B . Recently, Radler & Rheinhardt (2007) have 
pointed out that for finite values of B Q it may be important 
to consider instead the appropriate anisotropic expression, 
which can be written in the form 

TCijk (-u k Uj, p + p^bkbi 



The purpose of the present paper is to demonstrate that the 
values for both expressions, ([TJ and ©, are almost identical 
in the cases presented by BS05. We also show that the value 
of r, expressed in units of the turnover time, is in all cases 
close to unity, and in some cases better so than in BS05. 

2. Formalism 

As in BS05 we consider cases where the flow is driven ei- 
ther by a random body force in the momentum equation, or, 
alternatively, by random externally imposed currents in the 
induction equation. We calculated the isotropic expressions 



(3) 



in the presence of an imposed mean field Bo, where the tilde 
indicates the absence of the t factor, so a = t(<5k +<Sm)- As 
in BS05, we use additional superscripts k and m to indicate 
cases with kinetic or magnetic driving. The resulting values 
of ci^, a^, and presented below, are identical 

to those of BS05. In addition, we consider the appropriate 
component of the anisotropic expressions for the same sim- 
ulation data. Since in our case the mean field points in the y 
direction, and because we use volume averages and periodic 
boundary conditions in all three directions, we can write the 
anisotropic expressions for a yy in the form 
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where the superscript (a) indicates anisotropy. Again, we con- 
sider cases with kinetic and magnetic driving and thus obtain 
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Fig. 1. Dependence of <5^ ' and on R m in the kinetically Fig. 2. Same as Fig. [T] but for the relevant component of the 
forced case. Vertical bars give error estimates. (Adapted from 
BS05.) 
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and a^ 1 ^. The resulting 
values are normalized with respect to the corresponding rms 
turbulent velocities, 
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where u rms = [urmsMrms] 1 ^ 2 is the geometrical mean of the 
rms velocities for kinetically and magnetically driven runs. 
This particular normalization emerges naturally when deriv- 
ing the time scale r in Eq. (f2]). In the following we only con- 
sider the case of a statistically steady state, so b x b z 



u x u z , y , and hence also and a^', have converged to a 
stationary value. 

3. Results 

We consider the values of 5k and Sm and compare with the 
results of the appropriate component of the anisotropic ex- 
pressions; see FigsQ]and[2]for the kinetically driven case and 
Figs [3] and [4] for the magnetically driven case. The straight 
lines in Figs [T] and [3] denote fits to the data points, while 
in Figs |2] and [4] the same lines are just repeated as dashed 
lines and still represent only the fits to the isotropic data. This 
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helps demonstrating that the results change very little when 
the anisotropic expressions are used. 

It is remarkable that the differences between the isotropic 
and anisotropic expressions are rather systematic. Generally 
speaking, the anisotropic expressions give either the same or 
slightly smaller values than the isotropic expressions if the 
flow is driven hydrodynamically. The differences are larger 
for stronger fields (Bo = 0.1) and especially when the forc- 
ing it at larger scales (fcf=1.5). In that case the differences are 
around 15% and 25% for the kinetic and magnetic a effects, 
respectively. In the magnetically driven case the kinetic alpha 
effect tends to be smaller for the anisotropic expressions, but 
the magnetic alpha effect is either the same or larger for the 
anisotropic expressions. 

Following BS05, we also compare the results for all 
runs in tabular form; see Table Q] As in BS05, we non- 
dimensionalize the measurements for kinetically and mag- 
netically driven cases independently, because the root mean 
square velocities, tiJSs and i4ms, are different in the two 
cases; see Eq. 

There are two important aspects of the i? m dependence 
of kinetic and magnetic a effects. One is the fact that, at 
least for moderate values of R m , the two approach each other 
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Fig. 3. Dependence of and on R m in the magneti- Fig. 4. Same as Fig. [3] but for the relevant component of the 



cally forced case. Vertical bars give error estimates. (Adapted anisotropic expressions, a£ mJ and d^ mJ . The dashed lines 
from BS05.) represent the fit to the data of Fig. [5J not the present data! 



for finite field strength and increasing strength of the mean 
field. Furthermore, in the case of isotropic expressions, \&m\ 
could even slightly exceed the value of «k|- But when the 
anisotropic expressions are used, this is no longer the case-or 
at least less drastically so, e.g. in the middle panel of Fig. [2] 
The other aspect is the tendency for <5k to stay asymptotically 
independent of R m , even though the actual a effect decreases 
like with n = 0.5...1, as was shown in Fig. 2 of BS05 

for the same data. This property is critical to understanding 
the catastrophic quenching of the a effect for closed or peri- 
odic domains where magnetic helicity is a conserved quantity 
in the high conductivity limit. (We recall that, in contrast to 
the expressions for a$ and <3i/ , a itself was always calcu- 
lated as a = (£ ■ Bo) t /B , which does already account for 
the anisotropy for a. So the results for a remain unchanged 
from those obtained in BS05.) Let us also note in this con- 
nection that, within error bars, the off-diagonal components 
of the a tensor are found to be zero, i.e. | (£ x Bo)t | = 0. 

Finally we address the question of the relaxation time t. 
In BS05 we calculated r based on the values of a, a^ m \ and 

(k m) 

aiV . In the following we repeat the same analysis using the 
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. We recall that 



we allowed for different and unknown prefactors <?k an d 5m 
in front of oik and &u- We therefore wrote our unknowns in 
the form tc/k and t<7m> an d expressed them in normalized 
form as 

St#K,M = U Tms k{TgK,M- (6) 

These unknowns can be obtained by solving a matrix equa- 
tion which, in the present case, reads 

aWWaW 4 ak) Ust. 9K 

The result is shown in Fig. [5] for the old case using isotropic 
expressions of a, and in Fig.|6]for the present case using the 
anisotropic expressions. 

One of the most remarkable results from Fig. [6]is that the 
values of the magnetic and kinetic Strouhal numbers are in 
all three cases close to unity, whereas in the middle panel of 
Fig.[3]the Strouhal numbers were only about 0.3. In all other 
aspects the new results are rather similar to the old ones. For 
example, the values of magnetic and kinetic Strouhal num- 
bers are rather close to each other except in the case Bo = 0.1 
with k{ — 1.5, where the magnetic Strouhal numbers are 
somewhat larger than the kinetic ones. This is also the pa- 
rameter regime for which the largest differences were found 
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Table 1. Comparison of the results using the isotropic and anisotropic expressions for the various values of the normalized 
a for kinetically and magnetically forced runs. For fcf = 1.5 the resolution varies between 64 3 and 512 3 meshpoints for 
rj = 2 x 10~ 3 and 2 x 1CP 4 , corresponding to magnetic Reynolds numbers of 20 and 300, respectively, while for fcf = 5 the 
resolution varies between 32 3 and 256 3 meshpoints for r\ = 5 x 10 -3 and 5 x 10~ 4 , corresponding to magnetic Reynolds 
numbers of 4 and 60, respectively, The magnetic Prandtl number is always equal to unity, i.e. the viscosity v is always equal 
to the magnetic diffusivity, r\. 
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between Figs [T] and [2] Furthermore, like in BS05, we still 
find a drop in the Strouhal numbers in the case where R m is 
around 300. As argued in BS05, this may be connected with 
these simulations not having run for long enough. 

4. Discussion 

The work of BS05 was mainly an extension of earlier work on 
passive scale diffusion (Brandenburg et al. 2004), where cer- 
tain aspects of MTA were tested. In particular, it was shown 
that the relaxation time r in the r approximation is of the or- 
der of the turnover time (St = rtt rms fcf w 3). In the case with 
a magnetic field, the a effect was assumed to be expressible 
as a — t(<Sk + «m). The main result of BS05 was that St 
is independent of R m . This is important because neither (5k 
nor 5m decline with increasing values of R m . Instead, — (5m 
approaches <5k, resulting in near cancellation. Together with 
the finding that r is approximately independent of R m , this 
supports the validity of the assumed formula for a. It should 
be noted, however, that for R m ks 300 the result is not con- 
vincing and our present data suggest a drop in the Strouhal 
number. 

However, as RR07 have pointed out, several other issues 
remained open or unsatisfactory. In particular the compara- 
tive use of kinetically and magnetically forced models may 
be questionable. This was done to change the relative impor- 



tance of kinetic and magnetic a effects. The problem is that 
the nature of the turbulence can change considerably in the 
two cases. On the other hand, there is no reason why the ex- 
pressions for a should not apply equally well in both regimes 

Another problem is the use of isotropic expressions for 
(5k and <Sm- Surprisingly enough, as we have shown here, the 
isotropic expressions are indeed good proxies for the relevant 
component of the full anisotropic expressions. One advantage 
of using the anisotropic expressions is that the need for adopt- 
ing (slightly) different coefficients in front of (5k and (5m is 
now less severe, if at all present. 

Finally, there is the puzzle that, on the one hand, when 
using the first order smoothing approximation (FOSA), a is 
given by an expression involving just the actual velocity field 
while, on the other hand, according to the t approximation, 
it is the sum of magnetic and kinetic a effects. Obviously, 
a rigorous comparison between FOSA and r approximation 
is only permissible when the magnetic Reynolds number is 
below unity. In the present paper this is not the case, so the 
neglect of the higher order (triple) correlation terms under 
FOSA cannot be justified, given that the Strouhal numbers 
are always around unity. So this comparison may not have 
been permissible. However, the puzzle seems to exist even in 
the low magnetic Reynolds number limit, when the triple cor- 
relations can be neglected altogether. This case has been ana- 
lyzed recently by Sur et al. (2007), who showed that the for- 
mulations in terms of FOSA and t approximation are in fact 
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Fig. 5. Magnetic and kinetic Strouhal numbers as a function 
of R m for different values of Bq and kf. Here, kinetically and 
magnetically forced runs have been used to calculate sepa- 
rately gK ^ 5m- The horizontal lines are drawn to indicate 
the range over which the Strouhal numbers are approximately 
constant. (Adapted from BS05.) 
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Fig. 6. Same as Fig. [5] but for Strouhal numbers calculated 
from the expressions for the anisotropic alpha coefficients. 
The dashed lines represent the fits used in Fig. [5] and the 
solid lines represent new fits. 



equivalent (as they have to be, because the starting equations 
are the same!), but that the individual components contribut- 
ing to the total a-effect in the two formulations are different. 
In fact, it turns out that in the t approximation there is, in 
addition to the kinetic and magnetic alpha effects, in general 
also one more term resulting from the correlation between 
the small scale magnetic field and the forcing function. Only 
in the special case of ^-correlated forcing, that is adopted in 
many numerical investigations, does this extra term vanish. 
Nevertheless, even then the kinetic part of the alpha effect 
in the r approximation is not simply related to the alpha ef- 
fect obtained from the first order smoothing approximation, 
even if the actual velocity field is used in both cases. There- 
fore there is actually no puzzle in the limit of small magnetic 
Reynolds numbers either. 

5. Conclusions 

We have shown that the basic conclusions obtained in BS05 
carry over to the case where the anisotropic expressions for 



6iK and <5m are used. The present work provides an extra 
piece of evidence that the r approximation may provide a 
useable formalism for describing simulation data and for pre- 
dicting the behavior in situations that are not yet accessi- 
ble to direct simulations. There are currently no other ap- 
proaches capable of this. The basic hypothesis that the triple 
correlations are expressible as a damping term may not be 
well justified, although some important properties of this ap- 
proach seem to be borne out by simulations. A number of 
further practical tests of the r approximations could be en- 
visaged. One such example might be the so-called W x J 
effect of Rogachevskii & Kleeorin (2003, 2004), which was 
derived using the t approximation. Direct simulations of hy- 
dromagnetic turbulence with shear give qualitative support to 
this idea (Brandenburg 2005a), although it is not clear un- 
der which conditions the anticipated effect has the appropri- 
ate sign for dynamo action (Brandenburg 2005b; Riidiger & 
Kitchatinov 2006; Radler & Stepanov 2006). Further work in 
this direction would be worthwhile for establishing the real 
usefulness of the r approximation. 



511 



References 



References 



Acknowledgements. The Danish Center for Scientific Computing 
and the Center for Parallel Computers at KTH are acknowledged 
for providing computer time. 

References 

Bhattacharjee, A., Yuan, Y. 1995, ApJ, 449, 739 

Blackman, E. G„ Field, G. B. 2002, PRL, 89, 265007 

Blackman, E. G., Field, G. B. 2003, PhFl, 15, L73 

Brandenburg, A. 2005a, ApJ, 625, 539 

Brandenburg, A. 2005b, AN, 326, 787 

Brandenburg, A., Subramanian, K. 2005a, Phys. Rep., 417, 1 

Brandenburg, A., Subramanian, K. 2005b, A&A, 439, 835 

Brandenburg, A., Kapyla, P., Mohammed, A. 2004, PhFl, 16, 1020 

Field, G. B., Blackman, E. G. 2002, ApJ, 572, 685 

Gruzinov, A. V., Diamond, P. H. 1994, PRL, 72, 1651 

Kleeorin, N., Rogachevskii, I. 1999, PRE, 59, 6724 

Pouquet, A., Frisch, U., Leorat, J. 1976, JFM, 77, 321 

Radler,K.-H., Stepanov, R. 2006, PRE, 73, 056311 

Radler,K.-H., Rheinhardt, M. 2007, Geophys. Astrophys. Fluid Dy- 

nam. ( astro-ph/0606267 1 
Radler, K.-H., Kleeorin, N., Rogachevskii, I. 2003, GApFD, 97, 249 
Rogachevskii, I., Kleeorin, N. 2003, PRE, 68, 036301 
Rogachevskii, I., Kleeorin, N. 2004, PRE, 70, 046310 
Rudiger, G, Kitchatinov, L. L. 2006, AN, 327, 298 
Steenbeck, M., Krause, F, Radler, K.-H. 1966, Z. Naturforsch., 21a, 

369 

Sur, S., Subramanian, K, Brandenburg, A. 2007, MNRAS, 376, 
1238 



512 



